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Assuming that the order parameter corresponds to an equal spin triplet pairing symmetry state, we 
calculate the effect of phase fluctuations in quasi-one-dimensional superconductors at high magnetic 
fields applied along the y (b') axis. We show that phase fiuctuations can destroy the theoretically 
predicted triplet reentrant superconducting state, and that they are responsible for melting the 
magnetic field induced Josephson vortex lattice, above a magnetic field dependent melting temper- 
ature Tm- 



The upper critical field Hc^ of quasi-one-dimensional 
superconductors for perfectly aligned magnetic fields 
along the y [b') axis has been calculated (at the mean 
field level) for both singlet and triplet states In 
the case of triplet pairing, these studies predicted a 
remarkable reentrant superconducting phase, for mag- 
netic fields precisely aligned along the y (&') axis. New 
experiments performed in these systems lead to the 
observation of unusual superconductivity, for instance, 
in quasi-onc-dimcnsional superconductors (of the Bcch- 
gaard salts family with chemical formula (TMTSF)2X, 
where X = CIO4, PFg, ...) the experimental upper crit- 
ical field exceeds substantially the Pauli paramagnetic 
limit |4|j^. This indicates that high magnetic field su- 
perconductivity in these systems is most likely triplet. 
More recently, it has been argued by Lebed, Machida, 
and Ozaki |^ (LMO) that the spin-orbit coupling must 
be strong in order to explain the observed experimental 
upper critical fields j5| and the absence of the Knight 
shift for fields parallel to the y (6') axis. However, es- 
timates of the value of spin-orbit coupling are several or- 
ders of magnitude smaller than the values required to 
fit the critical temperature of (TMTSF)2PF6 §J|l, for in- 
stance, at low magnetic fields. Very recently, the present 
authors ||l0[| have analysed the angular dependence of 
iJc2 for magnetic fields applied along the yz plane, where 
weak spin-orbit coupling and equal spin triplet pairing 
were assumed. There, it was found that a strong sup- 
pression of the critical temperature Tc{H) occurred for 
very small angular deviations from the y (6') axis. This 
result is in qualitative agreement with the rapid drop of 
the putative upper critical in (TMTSF)2C104 |]. How- 
ever, all these theoretical works jl|-|^,|6|,|l0|| , have used a 
mean field approach, where the effects of fluctuations 
were completely ignored. Thus, the main point of the 
present paper is to discuss the effects of phase fluctu- 
ations on the field versus temperature phase diagram 
of quasi-one-dimensional superconductors at high mag- 
netic fields parallel to the y (b') axis. For that pur- 
pose, we study phase fluctuation effects in the weak 
spin-orbit, equal spin triplet state (ESTP) proposed by 
one us [pd|-p3[ as a possible candidate state for triplet 
superconductivity in the (TMTSF)2X family. In par- 
ticular, we show that phase fluctuations melt the pre- 



dicted magnetic field induced Josephson vortex lattice 
at high magnetic fields ||,|ll[l at a melting temperature 
Tm{H) < Tc{H). Furthermore, we show that the cur- 
vature of Tm{H) is opposite to that of Tc{H) at high 
magnetic fields (range from 5 to 20 T), and that the su- 
perconducting state still exists at high magnetic fields 
for temperatures T < Tm{H). 

We model quasi-one-dimensional systems via the en- 
ergy relation 

ea,o-(k) = eQ(k) - a^isH, (1) 

with the a-branch dispersion 

ea(k) = vpiakx — kp) + ty cos{kyb) + cos{kzC), (2) 

corresponding to an orthorombic crystal with lattice con- 
stants a,b and c [p^ along the x,y and z axis respec- 
tively. In addition, since Ep = vpkp ^ ty ^ tz, the 
Fermi surface of such systems is not simply connected, 
being open both in the xz plane and in the xy plane. 
Furthermore, the electronic motion can be classified as 
right-going (a = +) or left-going (a = — ). 

In this paper, we consider only magnetic field applied 
along the y (&') direction for simplicity, since quantum 
confinement along the transverse z direction occurs when 
tz/^c^ ^ 1, where ujc^ = vpGz with Gz — \e\HyC. 
Throughout the paper, we work in units where the 
Planck's constant and the speed of light are equal to 
one, i.e., Ti = c* = 1, and use the gauge A = [HyZ, 0, 0), 
where a and kx^ ky are still conserved quantities (good 
quantum numbers) while kz is not. 

In the presence of the magnetic field iJ, the non- 
interacting Hamiltonian is 

i7o(k- eA) = £„(k- eA) - cr/iijiJ (3) 

in the gauge A — [HyZ — Hzy, 0, 0). The eigenfunctions 
of i?o(k — eA) are 

'I>g„(r) = cxp[i{kxX + kyy)]JN,-n, (^) ' '•^^ 

where r — {x,y,z), z — rizC with associated quantum 
numbers qn = a,kx,ky, Nz, u and eigenvalues 

fgn = £a,a{^p) + CtNzUJc, ■ (5) 



The function Jp(u) is the Bessel function of integer order 
p and argument u, while now 

eQ,cr(kp) = vpiaK - kp) + ty cos{kyb) - (t^ibH (6) 

is a 2D dispersion. Notice that the electronic wave- 
function in Eq. is confined along the z direction 
when tz/oJc, ^ 1, thus limiting the electronic mo- 
tion to a nearly two dimensional situation. Further- 
more, notice that the eigenspectrum in Eq. involves 
many magnetic subbands labeled by quantum numbers 
and that the eigenvalue eq„ is invariant under the 
quantum number transformation (a, k^, ky, N^, cr) — + 
(—a, —kx, —ky, —Nz, cr), for the same spin state a. Since, 
these magnetic subbands are spin-split into spin-up and 
spin-down bands, Cooper pairs can be easily formed in a 
ESTP pairing state, involving electrons with quantum 
numbers {a,kx,ky, Nz,(7) and {—a,—kx,—ky,—Nz,a), 
provided that the pairing interaction conserves spin 
(which seems to be the case for the (TMTSF)2X family, 
except for very small spin-orbit and dipolar couplings) 
Thus, in the ESTP state the order parameter vector 



A(r) = [A+i(r),0,A_i(r)] 



(7) 



has components only in the ~ +1 {fi =tT) ^-^d mg — 
— 1 (/i =ii) channels, only. 

Now we turn our attention to the construction of the 
effective free energy, where we assume that the inter- 
action = A and the density of states TV^ = JV are 
independent of the spin channel, which implies that 
A^(r) = A(r). Following the functional integral for- 
mulation of the ESTP discussed by Sa de Melo |0, we 
derive the effective free energy functional in the limit 
where the magnetic field is exactly pointing along the y 
direction and where it nearly confines the electronic mo- 
tion to a two dimensional regime {tz/coc^ <C 1). In this 
case, A(r) = A„(x,j/), where n now labels the planes 
z = nc, and the effective free energy takes the form 

= H / (-^1 + -^2 + ^3 + ^4) + ^rn (8) 

The first term is the local quadratic form 

Ti^ai\Arf, (9) 

where the critical temperature can be derived from the 
coefficient 



ii = 



In 



T 



In 



(10) 



The second term has the form 

T2 = as. \dA,Jdxf + a2„ | {d/dy - i2Gyx/b) A„|2, (11) 



and corresponds to the spatial variation of A„(a;,y), 

4/ 



with coefficients = Af(3tGvp/T'^ and 02 



MC{tybf/2T'^, with C = 7C(3)/167r2. The third term 
corresponds to the contribution of the magnetic field in- 
duced Josephson coupling 



= a3|A„+iexp(-i2G2a;/c) - A, 



2 

n I I 



(12) 



with coefficient — J\f {tz/ujc^) In {\^c^ l7/27i'T). While 
the fourth term corresponds to the non-local fourth order 
contribution 



J^4 = a4i|A„|4 + a42|A„|2|A 



2 

n+l| 1 



(13) 



where 041 and 042 are complicated functions of ty, tz, 
ujr and ojr ■ The last term 

J Stt 

is just the magnetic energy. Using A„ = |A„| exp(i0„), 
it is easy to show that the terms in Eqs. (^ and ( p^ do 
not contribute to the phase fluctuation free energy. 



J~ p ^ ^ dv {^J~ px 4" J~ py ^" ^ pjs) ^" ^ prn- 



(14) 



Using units where h = c* = 1, the first two contributions 
to the phase fluctuation free energy are 



E„ 



i2eAr. 



(15) 



where the characteristic energies in the xy plane are 



= |a2j|A„|^ 
with fi = x,y. The third contribution is 



(16) 



^p,, = Jz COS {4),i+l 

(n+l)c 



where = 
term being 



„ - 2G.X - 2eA„ J . (17) 
dzAz/c, with the Josephson energy 



cV, = a3|A„||A 



n+l\ 



(18) 



The last contribution to the phase only fluctuation free 
energy is 



•^P-. = ^ / (V X Af , 



(19) 



where A is the fluctuation vector potential. The saddle 
point amplitude of the order parameter has the form 

|A„| = Ao [1 + itz/uc)^ cos(2G,a:)] , (20) 

where the prefactor is Aq = 5 [l + {tz/uJcJ'^] , (for Ep ^ 
wcj, with 6 = ^/Tc{H) [Tc{H) - T] /2C. The saddle 

point phase (l)n^ — 2Gznx — {^iTaj,^^! Hy)ns\n{2Gzx). 
This corresponds to a rectangular Josephson vortex lat- 
tice with periodicity 4 = i^/Gz {Gz is the inverse 
magnetic length) and Iz = c {c is the unit cell along 



the z direction) pj] ] and holds a flux quantum in- 
side the plaquette (Ixjz), i-e., i?44 — 4>o- In 
what follows we use similar methods to those developed 
by Horovitz (l^Jl^ and Korshunov and Larkin to 
study the vortex lattice melting in layered superconduc- 
tors. Writing down the phase of the order parameter as 

4>n{x, y) — 4>n^ + Xn{x,y), and integrating over the fluc- 
tuating vector potential we obtain the effective non-local 
sine-Gordon "Hamiltonian" 



H = J dr {Hi 



H2) 



/.dXn dXn 



(21) 



(22) 



(23) 



■f^2 = -r ^ COs(Xn-H - 2Xn + Xn-l), 

n 

where the co-sinusoidal couphng constant is 

r = j|/4/^^7,(7r)7j,(7r)/i2 (24) 



with the function 



7^(9) = 



A('Z) 



X{q) + 167T^E^d/^l 



(25) 



being the discrete Fourier transform of 7(71, ri'), where 
A(g) = [1 — cos(g)]. When T — 0, the effective Hamil- 
tonian given in Eq. (^ reduces to a layered anisotropic 
XY model, where phase fluctuations between layers n 
and n' are still coupled via the function 7^(n, n'). Scal- 
ing the integration variables x — > x^/E^jx and y 
y^Ey^iy, and defining the charges (j„, = /g^.. (ixn(r)/27r, 
which correspond to the vorticities at position r; in the 
n-th plane, leads to a partition function identical to the 
anisotropic quasi-two-dimensional Coulomb gas, when 
an expansion in powers of the very small parameter F 
and Gaussian integration over Xn are performed [p^ , p^ , 
i.e., 



Z = 



m i=0 \ iii / 



(- 

\2T 



2m 



exp(^), (26) 



where the exponent in the partition function is 



^ 2m 



(27) 



^,3 



where the charges qn 



1 for i — 1,....,TO and g,. 



— 1 for i ~ m + 1, 2m. In addition, = 
riij = Hj — rii, and the correlation function 



r, and 



M{r,j,nij) 



d^k dq ~ 



(27r)2 27r 



Af(k, q) exp [i(k • Vij + g)] 



has as Fourier transform the function Af (k, q) = 
T/ y/ExEyk'^j{q), with 7(9) = \/^xiq)lyiq)- Upon inte- 
gration over k and proper core regularization, the rele- 
vant part of the interaction term becomes 



M{rij,n,j) 



27Ty^E^y 



ln{rij/aa) 



(28) 



where a{nij) — J dq exp{inijq) /[2Tr'y{q)]. Notice that 
the charges interact logarithmically whether they are in 
the same layer or not; however, the strength of the log- 
arithmic interaction is larger when they are in the same 
layer. This means that in the limit where F ^ there 
is a phase transition which belongs to the same univer- 
sality class of the Kosterlitz-Thouless-Berezinskii tran- 
sition pc| , |2l| ], and corresponds, in the present problem, 
to the melting of the magnetic field induced rectangu- 
lar Josephson vortex lattice. The transition here occurs 
when 



T„ 



TT Ex{Tm)Ey{Tm) 



1 + ^Ex{Tm)Ey{T„,)16TT3c/(l)l 



(29) 



Notice that the right hand side of Eq. ^ is also de- 
pendent on temperature via E^{T). Using the expres- 
sions for E^ (T) defined in Eq. ( p^ and solving equation 
Eq. ( p9| ) for the melting temperature at infinite field re- 
sults in 



T',„(oo) 



TJ 00) 



1 



V 



(30) 



where rj = 2n^/2Tc{oo) /ty . Since 77 > 0, the melting tem- 
perature Tm(oo) is smaller than the saddle point (mean 
field) critical temperature Tc(oo), thus indicating that 
classical phase fluctuations reduce the transition temper- 
ature from Tj,(oo) to Tm{oo). However, this reduction is 
small for Bechgaard salts since rj 1. When the condi- 
tion — > 00 is relaxed the first order correction to 
can be computed using a perturbative renormalization 
group method ||l^,|l^,|l9| since E^ is much smaller than 
Tc{oo) and Tm{oo). This standard procedure leads to 



where e — 7r/2, — tt/uJc 
pendend coefficient 



2r„(oo) 



(31) 



and the magnetic field de- 



2ttTc{oo) 



Notice the explicit magnetic field dependence of F ~ 
H^^ \v?{H), such that the correction to rm(oo) increases 
with decreasing magnetic field as Tm{H) — Tm(oo) ~ 
H~^^\n'{H) in the high field regime where tz/ujc^ ^ 1. 
This behavior has the following interpretation: a reduc- 
tion of the magnetic field increases the magnetic field 



induced Josephson coupling and the system starts to 
become less two-dimensional, so that phase fluctuations 
become less efficient and Tm{H) increases with decreas- 
ing field. However, the increase in T„i{H) with decreas- 
ing field is very slow, and it looks quite flat when Tc{H) 
and Tm,{H) are plotted in the same scale (See Fig. |). 
This is in sharp contrast with what happens with the 
saddle point (mean field) critical temperature, which at 
high magnetic fields decreases from Tc (oo) as the mag- 
netic field is lowered. The difference in behavior between 
Tc{H) and T^{H) is illustrated in Fig. 0. 



that we have discussed here the effects of classical phase 
fluctuations for perfect alignment along the b' direction 
only. [H 
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FIG. 1. This figure shows the (mean field) critical tem- 
perature Tc{H) and the melting temperature Tm{H), in the 
regime tz/oJc^ <C 1 for applied fields H precisely aligned 
with the y {b') axis. The ratio t^/ujc^ « 1/3,1/18 for 
fields H = 5, 30 Teslas, respectively. Temperatures are in 
Kelvin, and fields are in Teslas. The parameters used were 
rc(oo) = 1.5K, ty — lOOK, tz — 5K with lattice constants 
characteristic of Bechgaard salts. 

In summary, we have assumed an ESTP state as 
a plausible candidate for triplet superconductivity in 
quasi-one-dimensional systems Q-^] , and we have pre- 
sented analytical results for the effects of phase fluctu- 
ations in the magnetic field versus temperature phase 
diagram of quasi-one-dimenional superconductors. We 
discussed for simplicity just the case of perfect align- 
ment with the y (b') direction, and showed that phase 
fluctuations destroy the reentrant superconducting phase 
when tz/iOc^ ^ 1- This loss of phase coherence cor- 
responds to the melting of the magnetic field induced 
Josephson vortex lattice above the melting temperature 
Tm{H). At very large fields Tm{H) is only slightly 
smaller than the mean field critical temperature Tc{H), 
however Tm{H) has the opposite curvature of Tc{H). 
This means that T„i{H) decreases with increasing field 
while Tc{H) increases with increasing field, provided 
that the condition tz/uJc^ <C 1 is satisfied. Furthermore, 
r,„(i?)-T,„(oo) ~ H-^Hii'{H) in this high field regime, 
and Tm{H) looks quite flat when plotted in the same 
scale as Tc{H). Finally, it is important to emphasize 
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